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Abstract-- Kelvin effect (Skin effect) is used in surface 
hardening produced by induction heating of gears, camforms, 
camshafts and other work pieces of fairly complex geometries. 
The induction heating equipment for surface hardening of metals 
and alloys (using LF or medium frequencies in the jargon of 
induction heating) is composed by a coil or coil assembly and a 
power semiconductor driving system up to 50kHz. The load seen 
by the driving system is equivalent to a transformer. The 
primary corresponds to the excitation coil or coil assembly, and 
the work piece correspond to a short-circuited secondary. To 
asses the electrical load it is necessary to determine the variations 
in skin depth from place to place due to local curvature effects in 
the work piece, and its variations in space and time due to 
variations in conductivity and magnetic properties coupled with 
thermal effects. In these and others technical applications of 
Kelvin effect it is often necessary to be able to relate local skin 
depths with local curvatures of the surface of electrically 
conductive bodies. The purpose of this paper is twofold. First, to 
derive a closed form analytical formula that relates the local skin 
depth with the local mean curvature and the well known skin 
depth for a flat conductive body. The limits of applicability of 
this formula are discussed. The predicted skin depths are 
compared with available experimental results obtained in the 
framework of surface hardening processes. Second, to apply the 
above mentioned formula to describe the electrical load of the 
induction heating equipment in the conditions used for surface 
hardening. In the choice or design of an induction heating system 
the parameters of the intended process (depth of Kelvin effect, 
temperatures to be reached and duration of the heating process, 
amongst others) put restrictions over the coils and the power 
driving system. To determine the best shape and size of induction 
coils or coil assemblies, the complex thermal and electromagnetic 
processes produced in the work pieces must be simulated using a 
computer code. The possible applications of the formula given in 
this paper, during the planning of these digital simulation is 
discussed. In the design of the electronic power system, driving 
the coils with the 
work pieces inside, the load impedance must be specified. The 
formula given here is applied to asses some characteristics of the 
load that may be of interest in the choice or design of the driving 
system. 


Index Terms-- Induction heating, Load modeling, Skin effect. 





I. INTRODUCTION 


Induction heating is a non-contact heating method. An 
electrically conductive work piece (a metal or a 
semiconductor) is located in an alternating magnetic field of a 
coil or coil assembly. By electromagnetic induction the field 
produces eddy — currents in the material, which is heated as a 
result of Joule effect. In ferrous metals, hysteresis loss heating 
occurs, but this effect is in general fairly small compared with 
the heating effect of the eddy- currents. The induced current 
density and its heating effects, are non uniform. The current 
density is maximum at the surface of the work piece nearest to 
the coil conductors. In the case of a simply connected body, 
the current density is zero at its centre. As the local density of 
thermal power is proportional to the square of the local 
electric current density, the heating effect is also non - 
uniform. If the frequency of the alternating magnetic field is 
high enough, the induced electric currents are developed in a 
thing layer adjacent to the above mentioned surface of the 
work piece. This is the skin or Kelvin effect, which is used in 
engineering manufacturing industry for the hardening of 
bearing surfaces, for welding, for brazing, and similar heat 
treatment processes [1]. 

The depth of penetration of eddy — currents can be 
characterized by the so called “skin depth”. It is a function of 
the frequency of the field, as well as, the physical properties 
and geometry of the work piece [1]. For surface hardening, the 
material is heated during a few seconds. In the earlier stages 
of the process, the heat is confined in a layer with the skin 
depth. The heated area is then cooled very fast, so that a 
surface hardening effect is obtained while the bulk of the 
material remains ductile. The distribution of heating in the 
body can be controlled adjusting the frequency of the 
alternating magnetic field. Frequencies in the range 50 Hz to 1 
MHz are used for several purposes, from through heating at 
low frequencies to surface heating at high frequencies. 

In this paper, surface hardening using power sources up to 50 
kHz will be considered. These are the so called medium 
frequency heating equipments. They are based on high — 
power thyristors. In some of these equipments, the rectified 50 
Hz supply is chopped by the thyristors and fed into a resonant 
circuit formed by the reactants of the loaded work coil 
(resistive and inductive) together with tuning capacitors. The 


load circuit seen by the driving system is equivalent to a short 
— circuited transformer. The primary corresponds to the 
excitation coil or coil assembly. The work piece can be 
considered as a short — circuited secondary. The impedance of 
this short circuited transformer depends on the geometry and 
physical properties of the work piece, the frequency of 
operation, and the geometry and number of turns of the coils. 
In the choice or design of an induction heating system to be 
applied in surface hardening processes, the parameters of the 
intended process (depth of Kelvin effect, temperatures to be 
reached and duration of the heating process, amongst others) 
put restrictions over the coils and the driving system. 

In surface hardening and other technical applications of 
Kelvin effect, it is often necessary to be able to relate the 
variations of the skin depth from place to place, due to local 
curvature effects in the heated body, and its variations in 
space and time due to variations in conductivity and magnetic 
properties coupled with thermal effects. A complete 
description of the problems posed by the method can be found 
in reference [2]. 

The practical importance of finding a relation between the 
local skin depth and the local curvature of the boundary of the 
work piece was stressed in reference [3], where an interesting 
qualitative description is given. 

The purpose of this paper is two fold. First, to derive a closed 
form analytical formula that relates the local skin depth with 
the local mean curvature and with the well known skin depth 
for a flat conductive body. Second, apply the above mentioned 
formula to describe the load impedance that must be driven by 
the induction heating equipment in the usual conditions for 
surface hardening of gears and other work pieces of fairly 
complex geometries. 


Il. AN ANALYTICAL FORMULA THAT RELATES SKIN DEPTH 
WITH LOCAL MEAN CURVATURE 


In order to study the Kelvin effect in low frequency fields, 
as in the case of induction heating, displacement current may 


be neglected [4]. In this case if S=EAH isthe Poynting 
vector, from Maxwell equations we obtain: 
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Let us consider the work piece located inside the working 
coil. If an alternating magnetic field is produced in the coil, 
the electromagnetic field in the air outside the work piece but 
inside the coil may be described by quasi - static equations. In 
the surface of the work piece the electromagnetic energy 
enters to the material along the normal to the considered point, 
due to the big difference in phase velocity between the air and 
the conductive body (refraction with absorption in the 
boundary) [5]. As a consequence, the Poynting vector in the 
boundary of the work piece is perpendicular to this surface 
and points towards the interior. 





V-S=-E-J-H- 





(1) 


> 


Let us consider now the lines of the field S inside the 
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work piece, and the congruence of surfaces orthogonal to 
these lines (one of these surfaces is the boundary of the heated 
body). But S= Sé, being E ç a unit vector and S the norm 


of the Poynting vector. One surface belonging to the above 


mentioned congruence pass through the point in which S is 
located, and é gis the normal vector to the surface at this 


point. According to a result of differential geometry and 
classical field theory applicable to a congruence of regular 
surfaces, the divergence of the normal vector is twice the 
mean curvature of the surface 49 [6] [7]. 


Then, 
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Here is the directional derivative of S in the 


direction of es , and s represents the arc length along the 
corresponding field line. From (1) and (2), substituting 
J = OF and taking time averages it follows: 
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where o is the conductivity of the material. If the magnetic 








(3) 





hysteresis is negligible, H DF may be put equal to zero. 
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Then it follows that: 
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ds 


The spatial scales that characterize the variation of the 
fields in a direction tangent to the boundary of the work piece 
are of the same order of magnitude of a characteristic 
dimension of the heated body. The spatial scale of variation of 
the fields along the direction of the Poynting vectors is of the 
same order of magnitude of the skin depth & for a flat 
boundary [4]: 

AE 
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Here @ is the angular frequency of the alternating magnetic 
field produced by the equipment and u is the magnetic 
permeability of the material. 

In the case of surface hardening, the skin depth ô is several 
orders of magnitude less than a characteristic dimension of the 
heated body. As a consequence, in a first approximation, the 
electric and magnetic fields can be described by the formulae 
corresponding to a flat boundary [4] [8]. Then 
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Here ( ) represents the RMS value of the magnitude, and 


1 
the factor V2 appears because the electric field has a phase 


lag of 7/4 relative to the magnetic field (neglecting magnetic 
hysteresis). Between the electric and magnetic fields there is 
the relation 
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where [4] [8] 
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From (4) to (7) it follows that: 
d (In S ) 2 
==2 8 
As 9 5 (8) 
If the boundary is flat, (9 = 0 and from (8): 
_28 
S(s)=S(O)e ê (9) 


gives the vanishing of the average magnitude of the Poynting 
vector going towards the interior of the body from its surface. 


S (0) is the magnitude of the average Poynting vector in a 


point of the surface of the work piece and S (s) is the 


corresponding magnitude at a point of the same field line at a 
distance s form the surface. This result suggests the 
introduction of an equivalent skin depth 6, that corresponds to 
the case of a curved surface: 
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This formula gives the relation between the local skin depth, 
the local mean curvature of the boundary and the skin depth 
corresponding to a flat boundary between the work piece and 
the surrounding air. As far as we know, this relation is a 
new one. 
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HI. A DESCRIPTION OF THE LOAD IMPEDANCE SEEN BY THE 
INDUCTION HEATING POWER SOURCE 


Now, let us consider the work piece as a short-circuited 
secondary of a transformer equivalent to the loaded coil. In 
the conditions prevailing during surface hardening the skin 
depth is much smaller than the geometric dimensions of the 
heated body. In this case we can consider a bar of non-circular 
cross section. The boundary of each cross section is given by 
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the same curve [ that may be described in the x-z plane of an 
orthogonal system of coordinates xyz. The y axis is parallel to 
the axis of the bar. Let h be the length of the bar. It is possible 
to introduce a surface density of electric current integrating 
the volumetric density in depth following the lines of the 


Poynting vector field. If J s is the phasor that represents the 


surface density, and Í is the global current that is induced in 
the work piece, then 


[= hd, (11) 


The electric currents flows parallel to the x-z plane, in each 
cross section of the bar. 


rE s is the phasor that represents the electric field tangent to 


the boundary of the heated body (also contained in the x-z 
plane), the induced electromotive force is given by: 


= pE,dl (12) 


In the case of a half infinite medium with a flat boundary, it is 
possible to show that [9] [10]: 
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If the Kelvin effect is strong enough, and the curvature 42 is 


(13) 


less than 1/6, we can suppose that (13) may be applied to the 
bar substituting 6 by the local skin depth 6, given by (10). 
Then multiplying and dividing by h and taking into account 
(11) it follows: 


z fÊ,al {a 
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So we obtain the following formula for the impedance of the 
work piece: 
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If IT) is the length of the curve T, from (10) and (15) we 


obtain 
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Here, by definition 
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This ô, is an average skin depth that corresponds to an 

average mean curvature: 
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=—— ped (18) 
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IV. DISCUSSION AND CONCLUSIONS 


In order to apply the results obtained above for the skin 
depth and for the impedance of a work piece, we will consider 
first the simplest and more common type of gear: spur gears 
[11]. Spur gears are used to transmit rotary motions between 
parallel shafts, with gears teeth that are located on the 
circumference of a circular plate. At the teeth level the 
boundary is convex. But if the boundary is convex, the mean 
curvature is negative. Then, the local skin depth, according to 
formula (10) will be greater than the skin depth of a flat 
boundary. If the boundary is concave, the mean curvature will 
be positive. Then the local skin depth will be smaller than the 
skin depth of a flat boundary. Now, let us suppose that this 
toothed wheel is located inside a coil in order to be heated by 
induction. The induced electric current will circulate in a layer 
adjacent to the surface of the teeth and the notches between 
the teeth, vanishing towards the interior of the wheel. As the 
local skin depth is higher at the teeth level than at the notches 
level, the real part of the local value of the impedance (the so 
called active resistance) will be lower at the teeth level than at 
the notches level. However, the same electric current 
circulates everywhere (so that the electric current density at 
the teeth level will be smaller than the electric current density 
at the notches level). As a consequence, the material adjacent 
to the notches will be heated more than the material of the 
teeth. This fact is important from the standpoint of induction 
hardening. Moreover, from the results derived in this paper it 
is possible to calculate the distribution of heat power in a layer 
adjacent to the surface of the work piece. This distribution 
may be used as an input in the equation of heat transfer , and 
solving this equation the corresponding temperature field may 
be calculated. 

If the frequency of operation is high enough, the skin 
depth for a flat boundary Ê will be so small that O will be 
negligible in comparison with 1. In that case, the local skin 
depth ó, will be equal to ô and the heat power produced by 


the induced electric current will be the same in all locations of 
the surface of the gear. 

These consequences of the analytical formulae are 
confirmed both by experimental results and by digital 
simulations of the process of surface hardening of gears by 
induction heating. [2] [3]. 

The assumptions made in order to derive (10) suggest that 
the formula for the local skin depth can be applied only if the 
product O”9 is small enough, perhaps less than 0.3. This 


guess may be substantiated in a comparison between the 
approximate analytical formulae obtained in this paper for 
bars of any cross sections and the exact results obtained for 
the Kelvin effect in infinite circular cylinders, in two different 
situations. In one of them, an infinite cylindrical conductor 
carries an alternating current parallel to its axis. The currents 
produce an alternating magnetic field perpendicular to the axis 
of the cylinder [3] [8]. In the other situation, an infinite 
cylindrical conductor is located in an alternating magnetic 
field, uniform and parallel to its axis. The induced currents are 
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contained in the circular cross sections of the cylinder. This 
case, studied by Förster and others [12], is a simplified 
mathematical model that may be applied in the framework of 
certain electromagnetic methods for non destructive testing of 
metals. However, it has a closer connection with the induction 
heating case. A detailed analysis of the relation between the 
approximate formulae obtained here and the consequences of 
the exact solutions for fields and currents in these cases will 
be done elsewhere. Here we only note that for a cylinder the 
mean curvature is constant and equal to twice the reciprocal of 


the radius 7} of the cylinder, so that from (10) it follows: 
ô 


=N (19) 
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Both exact solutions for the fields in the infinite cylinders, 
the classical one for the Kelvin effect in an alternating current 
conductor and the solution obtained by Forster for the 
conductor inside an infinite coil, are constructed using 
modified Bessel functions [8] [9] [12]. From an analysis of the 
asymptotic expressions of this Bessel functions, formula (19) 
is obtained, as will be shown in other work. 

Formula (16) gives an analytical estimate of the work piece 
impedance. If the line integral of the local mean curvature, 
taken over the curve’, is negative (as is the case of a 
predominantly convex body) both the resistive and the 
inductive parts of the impedance will be less than the 
corresponding impedance components for a flat body. If the 
line integral is positive (as is the case of a predominantly 
concave body) both the resistive and inductive parts will be 
greater than the corresponding impedance components of a 
flat body. If the frequency is high enough, these differences 
vanish. 

In order to calculate the impedance seeing by the electronic 
driving system, the work piece impedance must be reflected 
from the secondary to the primary circuit (the coil) of the 
equivalent transformer, using a well known procedure. End 
effects are neglected in this analytical approach, as well as the 
variations in conductivity and magnetic permeability produced 
by local temperature increase and magnetic hysteresis. 
Hysteresis may be described using a complex magnetic 
permeability in the formula for the skin depth [9]. One 
consequence is that the active resistance is now higher than 
the inductive reactance in the impedance given by 
(16).Magnetic saturation may be taken into account using a 
suitably smaller value for the magnetic permeability in the 
formulae for the skin depth and for the work piece impedance. 

Anyhow, to determine the best shape and size of induction 
coils or coil assemblies, in order to heat by induction the work 
pieces of complex geometry as found in practice, the 
analytical approach is clearly insufficient. The complex 
thermal and electromagnetic processes produced in the work 
pieces and the interaction between the work pieces and the 
coils must be simulated using computer codes. The used of 
these codes is not straight forward, and wrong results can be 


ô, = 


obtained from the numerical calculations. 

Nevertheless, the analytical approach affords general 
results that may be used as guide lines about what to expect 
and about what to search in relation with the results of the 
digital simulations. 

Also the analytical approach gives us useful lower and 
upper bounds for the values of parameters such as the local 
skin depths or the impedances that may be used to select 
induction heat equipment and to establish some main 
characteristics of the heating process. All this seems to 
deserve further study. 
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